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ABSTRACT 


Conventional  approaches  for  a  distributed  low  probabil¬ 
ity  of  detection  communications  system  with  a  large  number 
of  unique  transmitters,  and  a  single  or  a  few  receivers,  re¬ 
quire  receiver  complexity  proportional  to  the  number  of 
transmitters.  To  improve  efficiency  in  terms  of  receiver 
complexity,  two  alternative  designs  are  analyzed  and  com¬ 
pared  to  a  reference  receiver  whose  complexity  grows  line¬ 
arly  as  the  number  of  transmitters  increases.  The  first  al¬ 
ternative  system  groups  the  transmitters  into  clusters  whose 
pseudorandom  noise  codes  have  some  chips  in  common.  The  re¬ 
sulting  receiver  would  then  perform  two  stages  of  processes: 
identification  of  the  transmitting  cluster  and  received  bit 
detection.  The  total  number  of  processes  required  for  any 
given  transmitter  would  be  substantially  less  than  the  tra¬ 
ditional  receiver.  The  second  alternative  design  would 
utilize  a  common  long  spreading  code  and  a  shorter  cycli¬ 
cally  shifted  spreading  code  in  each  transmitter.  The  re¬ 
ceiver  utilizes  the  cyclic  shift  property  of  the  fast  Fou¬ 
rier  transform  to  recover  efficiently  both  the  identity  of 
active  receivers  and  the  data  sent  using  a  single  branch. 
The  complexity  of  the  two  proposed  systems  is  compared  to 
that  of  the  reference  system. 
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EXECUTIVE  SUMMARY 


Large  distributed  networks  often  have  more  transmitters 
than  receivers,  like  the  uplink  in  a  single  cellular  network 
cell.  Such  a  system  can  also  use  direct  sequence  spread 
spectrum  (DSSS)  in  order  to  achieve  a  low  probability  of  de¬ 
tection  (LPD)  and  multiple  access.  DSSS  systems  achieve  LPD 
performance  using  pseudo-random  (PN)  codes.  Such  coding  al¬ 
so  improves  anti  jamming  (AJ)  performance  and  mitigates  in¬ 
terference  between  transmitters. 

All  designs  considered  in  this  thesis  have  transmitters 
that  use  a  unique  PN  code  for  LPD,  AJ,  and  multiple  access. 
In  order  for  a  receiver  to  process  all  signals,  it  must  have 
memory  to  store  the  reference  PN  code(s),  a  mixer  to  de¬ 
spread  the  signal,  a  matched  filter  to  extract  decision  da¬ 
ta,  a  sampler  to  select  discrete  data  for  analysis,  and  a 
decision  algorithm  in  hardware  or  software  to  make  bit  deci¬ 
sions.  In  a  traditional  receiver,  this  requirement  results 
in  significant  growth  in  complexity  as  the  number  of  trans¬ 
mitters  increases,  making  large  systems  unwieldy  to  field 
and  difficult  to  scale  for  increased  traffic  loads. 

One  potential  solution  involves  sub-dividing  the  re¬ 
ceived  signal  set  into  smaller  groups  using  a  mask  of  common 
PN  code  points.  For  example,  a  system  with  32  transmitters 
requires  32  copies  of  the  receiver  processing  hardware,  or  a 
software  radio  equivalent,  including  a  demodulation  algo¬ 
rithm  run  32  times.  Dividing  the  set  into  four  subsets  of 
eight  related  signals  improves  the  efficiency  of  processing, 
as  a  single  signal  requires  only  a  portion  of  the  receivers 
for  demodulation.  The  receiver  would  first  determine  which 

xiii 


subset  contains  the  received  signal,  then  attempt  to  demodu¬ 
late  all  signals  in  that  subset.  The  total  number  of  proc¬ 
esses  would  be  12  vice  32— four  initial  detection  decisions 
to  determine  the  subset  and  eight  demodulations.  This  is 
potentially  efficient  if  the  number  of  active  transmitters 
is  a  small  fraction  of  the  total  transmitters  at  any  given 
moment.  However,  in  the  case  of  all  transmitters  sending 
nearly  simultaneously,  the  coding  mask  receiver  is  less  ef¬ 
ficient  as  32  demodulations  and  four  detections  require  36 
processes,  while  the  linear-growth  receiver  still  only  re¬ 
quires  32.  Additionally,  the  probability  of  error  in  the 
coding  mask  system  will  be  worse  than  for  the  reference  re¬ 
ceiver  as  the  mask  detection  stage  introduces  additional  er¬ 
ror  . 

A  second  proposed  solution  involves  using  short  PN 
codes  related  to  a  single  code  by  cyclic  shift  for  multiple 
access.  To  provide  for  code  division  multiple  access 
(CDMA) ,  each  transmitter  is  preprogrammed  with  a  predeter¬ 
mined  shifted  form  of  a  base  PN  code  with  duration  of  a  sin¬ 
gle  bit.  Every  transmitter  then  uses  a  common  longer  PN 
code  to  provide  transmission  security  by  mixing  the  two 
codes  together  prior  to  mixing  with  data.  The  longer  PN 
code  has  duration  greater  than  a  single  bit  and  potentially 
never  repeats  within  the  duration  of  the  entire  message  the¬ 
reby  ensuring  a  unique  PN  code  for  each  transmitter  and  LPD. 

This  proposed  receiver  also  removes  the  common  PN  code 
through  mixing  before  filtering  and  sampling  of  the  combined 
received  signal.  It  compares  the  discrete  Fourier  transform 
(DFT)  processed  signal  received  against  a  reference  DFT  of 
the  short  unshifted  PN  code,  extracts  the  transmitter  iden- 


xiv 


tity  through  an  inverse  DFT  (IDFT)  and  completes  bit  deci¬ 
sions.  This  process  achieves  the  same  functionality  as  the 
matched  filter  in  the  reference  receiver,  but  with  less  com¬ 
ponents  required. 

The  IDFT  output  identifies  which  transmitters  sent 
data,  as  well  what  data  was  sent.  Increasing  the  number  of 
transmitters  requires  additional  processing  of  the  IDFT  out¬ 
put,  but  this  growth  is  smaller  as  the  system  needs  fewer 
copies  of  shift  size  and  bit  detection  components.  Finally, 
the  cyclic  PN  code  system' s  improved  scalability  is  not  at 
the  expense  of  error  performance.  It  has  equivalent  error 
performance  to  the  reference  receiver  for  binary  phase  shift 
keying  (BPSK)  modulation  within  an  additive  white  Gaussian 
noise  (AWGN)  channel. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

This  thesis  examines  potential  solutions  to  the  problem 
of  scalability  for  a  widely  distributed  direct  sequence 
spread  spectrum  (DSSS)  communications  system.  This  system 
must  have  low  probability  of  detection  (LPD)  characteris¬ 
tics,  which  prevent  the  exploitation  of  its  signals  by  an 
interceptor,  as  well  as  prevent  interference  between  trans¬ 
mitters.  Each  transmitter,  therefore,  must  have  a  unique 
pseudo-random  noise  (PN)  code,  and  that  in  turn  requires  the 
receiver  to  process  all  possible  system  PN  codes  simultane¬ 
ously. 

Adding  extra  copies  of  the  processing  hardware  for  each 
additional  transmitter  quickly  leads  to  an  unwieldy  system. 
There  will  be  growing  hardware  costs  and  physical  space  re¬ 
quirements,  as  well  as  potential  increased  delays  in  the 
signal  processing.  Throughout  this  thesis,  the  system  that 
requires  an  additional  set  of  hardware  for  each  additional 
transmitter  will  be  referred  to  as  the  "reference  receiver" 
and  its  detailed  analysis  is  contained  in  Chapter  III.  The 
reference  receiver  will  provide  the  error  performance  bench¬ 
mark  for  comparison  of  scalability  solutions.  The  reference 
receiver  and  all  solutions  to  the  scalability  problem  will 
use  DSSS  with  unique  PN  codes  in  each  transmitter  to  main¬ 
tain  LPD  and  multiple  access.  Additionally,  each  analysis 
will  use  an  additive  white  Gaussian  noise  (AWGN)  channel 
with  no  forward  error  correction  coding. 
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Low  probability  of  interception  (LPI)  and  low  probabil¬ 
ity  of  exploitation  (LPE)  are  sometimes  used  interchangeably 
with  LPD.  LPD  in  the  context  of  this  thesis  refers  to 
structuring  the  signal  such  that  it  is  relatively  immune 
from  being  detected  [1],  and  this  term  will  be  used  from  now 
on  in  discussing  preventing  the  exploitation  of  the  signal. 

B.  OBJECTIVE 

This  thesis  will  examine  and  validate  the  reference  re¬ 
ceiver,  as  well  as  analyze  two  potential  solutions  to  the 
scalability  problem  by  comparing  their  error  performance  and 
scalability.  One  solution  is  the  "cyclic  PN  code  receiver." 
It  uses  two  PN  codes  to  provide  multiple  access  and  identi¬ 
fication  of  transmitters,  as  well  as  LPD.  A  common  receiver 
with  a  single  branch  of  hardware  processes  all  signals  si¬ 
multaneously  using  discrete  Fourier  transforms  (DFTs) ,  which 
act  as  a  matched  filter  with  a  single  bit  decision  block. 

An  alternative  solution  is  the  "coding  mask  receiver." 
It  uses  the  same  hardware  as  the  reference  receiver  for  bit 
decisions  but  sub-divides  them  into  groups  with  similar  PN 
codes.  Each  subset  then  has  a  detection  branch  to  reduce 
the  total  number  of  processes  required  for  a  single  received 
signal.  This  solution  speeds  up  processing  of  communica¬ 
tions  by  requiring  fewer  processes  for  a  single  signal  re¬ 
ceived.  However,  the  same  physical  hardware  is  required  to 
support  the  full  system,  so  scalability  is  still  potentially 
an  issue. 
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c. 


RELATED  WORK 


PN  codes  with  different  lengths  in  conjunction  with  the 
discrete  Fourier  transform  provide  a  method  to  efficiently 
emulate  a  set  of  matched  filters  and  maintain  bit  error  rate 
performance.  The  cyclic  PN  code  receiver  solution  is  simi¬ 
lar  to  cyclic  code  shift  keying  (CCSK) ,  a  technique  used  in 
the  Joint  Tactical  Information  Distribution  System  (JTIDS) 
[2],  except  that  CCSK  uses  the  cyclic  shift  to  indicate  a 
data  value  while  in  this  thesis,  the  cyclic  shift  provides 
multiple  access  and  identifies  the  source  transmitter  at  the 
receiver  [3] .  The  coding  mask  receiver  solution  uses  unique 
PN  codes  with  points  of  similarity  to  perform  the  initial 
signal  detection  and  then  separate  subgroups  of  receivers 
for  demodulation.  In  terms  of  the  signals,  this  is  the  same 
as  a  typical  CDMA  [4],  but  with  prescribed  correlation  be¬ 
tween  portions  of  the  spreading  PN  codes. 

Additionally,  the  cyclic  PN  code  receiver''  s  shorter  PN 
code  repeats  at  the  bit  rate  and  provides  code  division  mul¬ 
tiple  access  (CDMA)  similar  to  the  methods  described  by  Lee 
and  Miller  and  separately  by  Ha  as  used  in  the  IS-95B  cellu¬ 
lar  standard.  There,  the  long  PN  code  scrambles  the  data 
while  two  short  codes,  with  different  shifts,  maintain  mul¬ 
tiple-access  within  and  between  cells  while  minimizing  in¬ 
terference  [4],  [5]. 

The  next  chapter  provides  an  overview  of  the  design 
points  of  the  various  receivers  investigated  in  this  thesis. 
Chapter  III  then  provides  a  detailed  analysis  and  perform¬ 
ance  simulation  of  the  reference  receiver.  Chapter  IV  fol¬ 
lows  with  a  detailed  description,  theoretical  calculations, 
and  simulation  of  the  cyclic-PN  code  receiver,  followed  in 
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Chapter  V  by  a  similar  analysis  of  the  coding  mask  receiver. 
Chapter  VI  provides  a  summary  of  all  data  and  compares  error 


rate  and 
provided 


scalability  performance.  Finally,  conclusions  are 
in  Chapter  VI I . 
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II.  SYSTEM  DESIGN 


In  order  to  operate  as  a  distributed  communications 
system,  a  proposed  design  must  meet  three  requirements.  It 
must  maintain  low  probability  of  detection  (LPD) .  This 
means  preventing  the  detection  of  the  signal  by  an  intercep¬ 
tor  through  the  design  of  the  signal's  waveform  and  does  not 
refer  to  protecting  the  data  through  encryption.  Next,  a 
proposed  design  must  have  error  performance  similar  to  un¬ 
coded  BPSK  in  an  AWGN  channel.  Superior  scalability 
achieved  through  an  increase  in  the  bit  error  rate  is  not 
acceptable.  Finally,  each  valid  solution  must  have  superior 
scalability  to  the  reference  receiver,  which  is  expected  to 
have  linear  growth  as  the  number  of  transmitters  increases. 

A.  LOW  PROBABILITY  OF  EXPLOITATION  DESIGN 

Direct  sequence  spread  spectrum  (DSSS)  addresses  issues 
of  LPD  by  spreading  a  signal's  spectrum  to  decrease  the 
probability  of  detection  [1].  Additionally,  using  multiple 
pseudo-random  noise  (PN)  codes  allows  for  increased  protec¬ 
tion  as  well  as  improvements  in  demodulation  as  each  re¬ 
ceiver  can  use  unique  inner  codes. 

1.  Non-recursive  Coding 

Effective  PN  codes  must  have  a  random  appearance.  Each 
transmitter  must  therefore  use  a  non-recursive  code,  i.e., 
one  that  does  not  repeat  for  the  duration  of  each  transmis¬ 
sion  and  ideally  never  use  the  same  code  in  subsequent 
transmissions.  Mixing  two  PN  codes  achieves  this  as  long  as 
one  of  the  codes  is  non-recursive  [8] . 
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A  random  number  generator,  as  simulated  in  MATLAB,  can 


generate  an  infinite  length  code  for  this  purpose.  In  a  si¬ 
mulation,  a  common  code  generator  can  provide  codes  for  the 
transmitter  and  receiver.  In  a  real  system,  synchronizing 
the  generators  with  a  common  seed  achieves  the  same  [9], 
[10]  . 


2 .  Low  Power 

An  LPE  communications  system  typically  uses  less  power 
in  its  transmitters  in  order  to  reduce  the  chance  of  detec¬ 
tion.  However,  without  forward  error  correction  coding, 
this  will  also  reduce  error  performance  and  range.  For  the 
models  in  this  thesis,  the  energy  per  bit  remains  constant 
while  the  noise  power  varies,  allowing  for  the  testing  of  a 
range  of  signal  to  noise  ratios. 

3 .  High  Bandwidth 

Spread  spectrum  signals  require  more  bandwidth  than 
non-spread  signals  to  transmit  at  the  same  data  rate  [5]  . 
The  proposed  solutions  control  mutual  interference  through 
PN  codes.  Future  design  work  should  explore  forward  error 
control  coding  to  support  binary  phase  shift  (BPSK)  or  quad¬ 
rature  phase  shift  keying  (QPSK) .  This  thesis'  simulations 
used  BPSK  modulation  for  simplicity.  However,  DSSS  designs 
using  QPSK  have  superior  LPE  performance,  since  the  signal 
squared  does  not  reveal  the  carrier  frequency  [1] . 

B.  SCALABILITY 

Examining  a  system  with  N  transmitters,  for  a  linear- 
growth  receiver  to  processes  the  received  signals,  it  must 
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repeat  the  filtering  and  bit  decisions  for  all  signals  N 
times,  resulting  in  an  approximate  total  number  of  multipli¬ 
cation  processes  of  N2  .  However,  the  fast  Fourier  trans¬ 
form  (FFT)  algorithm,  a  means  to  implement  the  discrete  Fou¬ 
rier  Transform  (DFT)  ,  requires  0.5Nlog2N  complex  multiplica¬ 
tion  operations  [6]  .  The  proposed  receiver  will  have  one 
FFT  and  one  inverse  FFT  (IFFT)  each,  requiring  a  potential 
total  of  N\og.,N  complex  multiplications,  which  increases  at 
a  rate  less  than  the  reference  receiver  as  the  system  grows 
in  size.  A  detailed  scalability  discussion  is  included 
Chapter  V,  section  B.  The  next  chapter  describes  the  refer¬ 
ence  receiver  in  detail. 
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Ill .  REFERENCE  RECEIVER 


The  reference  receiver  uses  a  unique  PN  code  for  each 
transmitter.  Its  receiver  has  some  common  components  that 
process  all  received  signals  together,  but  it  relies  on  sep¬ 
arate  de-spreading,  filtering,  and  bit  decision  blocks  for 
each  transmitter.  The  reference  receiver  is  representative 
of  a  straightforward  design  approach,  and  therefore  provides 
a  good  benchmark  for  comparisons.  Determining  the  reference 
receiver' s  BPSK  error  performance  and  validating  theoretical 
calculations  with  simulation  are  the  basis  for  comparison 
and  analysis  of  the  two  alternative  designs.  The  block  dia¬ 
gram  for  this  reference  design  was  based  on  Nicholson' s  hy¬ 
pothetical  DSSS  transmitter  with  a  single  PN  code  and  car¬ 
rier  frequency  oscillator  [1] . 

A.  REFERENCE  RECEIVER  DESIGN 


The  reference  transmitter  and  receiver  use  a  single 
spreading  code  and  upconversion  to  create  the  spread  spec¬ 
trum  waveform.  The  receiver  complexity  grows  linearly  as 
the  number  of  transmitters  increases,  but  it  provides  the 
benchmark  for  error  performance. 


Its  probability  of  bit  error  is 


0 


f 


V 


where  Eh  is 


energy  per  bit  and 


is  the  two-sided  noise  power  spec¬ 


tral  density,  which  is  identical  to  a  binary  phase  shift 
keyed  (BPSK)  system  [7]  and  derived  in  equations  (3.16) 
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and  (3.17)  .  This  uncoded  BPSK  error  performance  is  as  ex¬ 
pected  for  a  DSSS  system  because  the  signaling  is  antipodal 
[8]  . 

1 .  The  Reference  Transmitter 

The  reference  transmitter  initially  transforms  the  po¬ 
lar  binary  data  stream  jt/J  into  an  antipodal  non-return-to- 

zero  (NRZ)  binary  signal  mp{t)  at  the  bit  rate  Rh  = 

where  Pn(t-iTb)  is  a  pulse  function  for  the  i,h  bit  time  with 
duration  Tb  ,  and  subscript  p  in  mp  (?)  is  the  transmitter 
number  as  shown  in  equation  (3.1)  . 

oo 

mp(t)  =  Yjd:Prb(t-iTb)  (3.D 

-00 

The  data  mixes  with  the  transmitter' s  unique  pseudoran¬ 
dom  noise  (PN)  code,  cp  (7) ,  and  carrier  wave,  cos(2^/f^)  . 
Each  transmitter  uses  the  same  carrier  frequency  ( fc ) ,  which 
is  much  higher  than  the  bit  rate.  The  amplifier  applies  a 
preset  gain  of  A  before  the  signal  enters  the  channel  as 
s(7)  =  Am  (t)c  [t) cos [2 7rfct) ,  as  shown  in  Figure  1. 
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2. 


The  Reference  Receiver 


The  channel  for  this  model  is  assumed  to  be  additive 
white  Gaussian  noise  (AWGN) ,  adding  a  white  Gaussian  noise 
factor  n{t )  to  the  combined  signal,  resulting  in  a  received 

signal  of  s(V)  +  m(7)  as  shown  in  equation  (3.2)  .  At  the  out¬ 
put  of  the  LNA,  the  signal  received  still  matches  equation 
(3.2)  for  theoretical  calculations. 

r{t )  =  Amp  (t)cp  (7)  cos  (2^/^)  +  n{t)  (3.2) 

One  complete  receiver  branch  is  required  for  each  of  N 
transmitters,  resulting  in  a  system  as  shown  in  Figure  2. 


Figure  2 .  Reference  Receiver  System 


Examining  just  the  signal  received  from  the  pth  trans¬ 
mitter,  the  received  signal  passes  through  a  low-noise  am¬ 
plifier  (LNA)  before  mixing  with  the  carrier  and  unique  PN 
code  matching  its  transmitter.  The  signal  now  has  two  com¬ 
ponents,  the  squared  cosine  function,  and  noise. 
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(3.3) 


r{t)  =  Amp  (t)cp(t)cp(t)  cos  (2  xfct)  cos  (2xfct)  +  "{t)cp  (t )  cos  {2nfct) 
=  Amp  (i t)cp{t ) cos2  (2 nfct)  +  n(t)cp(t) cos (2 nfct) 


Using 

the  trigonometric 

identity 

c°s2(x)  =  ^( 

l-cos(2x)) 

[11],  and 

the  squared  PN  code 

equating 

to  one. 

equation 

(3.3)  simplifies,  giving  a  signal  with  three  components:  a 
signal  term,  a  double  frequency  term,  and  a  noise  term. 

a{t)  =  jmp(t)(l-cos(4xfct))  +  n(t)cn(t)cos(27rfct) 

2  (3.4) 

A  A 

=  ^mp(t)-^mp(t)COS(4xfct)  +  n(t)Cn(t)COS(2xfct) 

Looking  at  only  the  signal  term  entering  the  matched 
filter,  its  output  is 

A  A  °r 

y(t)  =  h(t)*  —  mp  (7)  =  —  J  pT(ya^mp(t  -  a^da  .  (3.5) 

^  -oo 

With  substitution,  the  time  interval  of  the  convolution 
integral  changes  from  -oo  <  a  <  oo  to  0  <a<Tb.  Using  the  new 

interval,  the  pulse  function  term,  ,  in  the  convolution 

integral  equates  to  one,  giving  a  new  form  of  the  integral. 

T 

A  oo  1b 

y(t)  =  —  ^  ;  J  \pT  (f  -  a  -  iTh  yia  (3.6) 

2-«,  o 

Evaluating  this  integral  at  the  sample  instant,  t  =  kTb , 

where  ke  Z,  yields  the  signal  portion  of  the  decision  sta¬ 
tistic. 

A  00 

y{kTb)  =  i:Yjdi\  PT{{k-i)Tb-a)da  (3.7) 

^  -oo  o 

In  evaluating  the  convolution  integral  with  k,ie  Z, 
there  are  three  possible  relations:  k—i> 2,  k—i< 0,  and 
k-i  =  1.  In  the  first  instance,  the  integral  is  zero  in 
0  <a<Tb.  The  second  condition  likewise  solves  to  zero  for 
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the  same  reason.  In  the  third  instance,  evaluating  the  in¬ 
tegral  at  k—i= 1  results  in  the  final  form  of  the  signal  of 
interest  at  the  output  of  the  matched  filter.  The  final 
value  shown  in  equation  (3.8)  is  the  mean  for  the  decision 
statistic  normal  random  variable  after  sampling  and  prior  to 
the  bit  decision. 


y(kTb) 


=  jX' dA-i.ij  Pt  ( Th-a)da 


A  r 
=  ~dk- ij  pT(Th-a)da 


A,  _ 

—  2  dk-\Tb- 


(3.8) 


3 .  Double  Frequency  Term 


The  double  frequency  term,  —  mp  (7)  cos  (4^/^) ,  convolves 

with  the  matched  filter  impulse  response,  h(t)  =  pT(t),  to 
yield 


A 


(3.9) 


=  —  J  pT[a)m(t -a)cos(47tfc(t -a)}da. 


Using  the  same  substitution  for  the  data  signal  as  in 
the  signal  of  interest  above,  the  interval  changes  from 
-oo<a<oo  to  0  <a<Th,  giving  a  new  convolution  integral  of 
the  form: 

T 

A  00  b 

y(V)  =  — J pT(t-a-iT)cos(yAnfc(t-ay)da  .  (3.10) 


Evaluated  at  t  =  kTh ,  where  keZ,  the  integral  becomes: 


A  00  b 

( kTb )  =  -X^i  Pr{(k-i)Tb  ~ a) cos (4xfc  (kTb  -a))da  . 


(3.11) 


0 
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Similar  to  the  signal  of  interest,  the  double  frequency 
convolution  integral  only  produces  non-zero  output  when 
k—i= 1,  resulting  in  the  final  form  shown  in  equation 
(3.12)  . 

A  Tr 

y{kTb )  =  -zdk-i  J  Pt  (Tb  -  a)  cos  (4  nfc  ( kTb  -a))  da 

^  0 

A  Th 

=  —dk  , | cos (4 nfc[kTh  -a)^da  (3.12) 

2  o 

=  77[8in(4« )  -  sin  (4*/e  ((*  - 1)  z; )) 

°^Jc 

The  term  in  square  brackets  in  the  solution  has  a  max 
value  of  +2  and  minimum  of  -2.  Therefore,  the  absolute 
value  of  the  solution  must  be: 

^^[sm(4^fckTb)-sm{47rfc((k-l)Tb))j  •  (3.13) 

Accounting  for  fc»\  and  its  location  in  the  denomina¬ 
tor,  all  possible  results  of  this  term  are  small  compared  to 
the  possible  results  of  the  data  signal  in  equation  (3.8), 
so  the  double  frequency  term  is  ignored  as  insignificant  in 
further  analysis  of  the  bit  decision  statistic. 

4 .  Noise  Calculations 


The  noise  term  enters  the  match  filter  as 
n(t)cn  (7)cos(2;r/d )  .  However,  analyzing  the  resulting  power  of 
the  noise  is  more  useful.  The  power  spectral  density  (PSD) 


of  noise  from  the  AWGN  channel  at  the  antenna  is  S  (f)  =  — -. 

nn  \  J  )  ^ 

Mixing  with  the  PN  code  has  no  effect  on  the  resulting  PSD. 
The  instantaneous  power  of  a  cosine  function  is  cos2(x). 
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which  has 
rier  wave 
the  power 


an  averaged 
results  in 
density  or 


of  — .  Mixing  the  noise  with  the  car- 
2 

a  noise  process  with  a  PSD  with  half 
,  x  N 

Sn„(f)  =  —^-r  which  is  the  PSD  at  the  in¬ 


put  of  the  matched  filter.  The  output  of  the  matched  filter 
has  power  spectral  density 

5.,(/)  =  5„.(/)|ff(/)|2=^|ff(/)|2  ■  (3.14) 

To  integrate  equation  (3.14),  Parseval ' s  theorem  [12] 
allows  the  time-domain  impulse  response  function,  h(t)  =  pT{t), 

to  be  used  instead.  Therefore,  the  power  of  the  noise  at 
the  output  of  the  matched  filter  is 


-oo 

(3.15) 

IMO  \dt 


Since  pT(t)  =  \  for  0  <t<Ths  and  has  zero  value  otherwise, 

equation  (3.15)  resolves  to  the  noise  power  at  matched  fil¬ 
ter  output. 


5 .  Reference  Receiver  Error  Performance 


Using  the  output  of  the  matched  filter,  the  signal  of 
interest  and  noise  enter  the  decision  block  where  the  re¬ 
ceiver  makes  bit  determinations.  Using  the  signal  of  inter¬ 


est  input. 


y(kTb) 


A 

2  Clk-lIb  ' 


the  decision  statistic  is  a  Gaus- 
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f  Ad  j1  NT) 

sian  random  variable,  N  — k~l  b , — —  1.  This  is  used  to  de- 

l  2  4  J 

termine  the  probability  that  the  decision  block  output  bit, 
dk  ,  is  in  error  when  the  sent  bit  is  dk_x=—  1.  Dividing  by 


the  standard  deviation  of  the  decision  statistic. 


one  finds  the  area  under  the  right-hand  tail  of  the  Gaussian 
distribution,  resulting  in 


Pr(<4- 1  = 1 1  dk- 1  =  -1)  =  p  *(y{kTb)  >  o) 


(3.16) 


Using  the  equivalence  A2Tb=2Eh,  where  Eh  is  the  energy 

per  bit,  the  resulting  probability  of  bit  error  is  identical 
to  BPSK,  as  expected  since  the  signal  is  antipodal  [7] . 


^{ddec=l\dk_x=-\)  =  Q 


2F 

=  Q  — - 

\  -AT 


(3.17) 


The  theoretical  error  performance  for  an  antipodal  sys¬ 
tem  plotted  using  a  range  of  signal  to  noise  (SNR)  ratios  or 
E. 

— —  (bit  energy  to  noise  power  spectral  density)  ranging 
N0 

from  0  dB  to  10  dB  is  shown  in  Figure  3.  This  curve  repre¬ 
sents  the  theoretical  best  error  performance  of  the  refer¬ 
ence  receiver  with  uncoded  BPSK  modulation  and  was  generated 
using  MathCAD  (v. 14 . 0 . 0 . 163) . 


1  The  notation 


N(X ,&2) 


is  used  in  this  thesis  to  describe  a  random 


variable  as  a  normal,  i.e.,  Gaussian,  random  variable  with  mean  X  and 
variance  <J  . 
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Figure  3.  Optimum  Bit  Error  Rate  for  a  BPSK  System  over  a 

E, 

Range  of  — — 

N0 

B.  MODELING  THE  REFERENCE  RECEIVER 

The  reference  receiver  was  simulated  with  MATLAB  and 
Simulink.  The  MATLAB  file,  shown  in  Appendix  A,  used  the 
variable  "sim_length"  to  control  number  of  bits  simulate. 
By  setting  this  variable  to  three,  the  simulation  ran  for 
100,000  bits  for  each  SNR  level  over  the  range  of  0  dB  to  9 
dB .  The  AWGN  channel's  calculations  were  independent  for 
each  bit. 
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1. 


Simulation  Results 


The  reference  receiver' s  simulation  results  are  shown 
in  Figure  4.  The  black  line  represents  the  BPSK  theoretical 
bit  error  rate  performance,  and  is  the  same  curve  from 
Figure  3.  The  simulation  bit  error  rate  results  are  the  red 
circles  on  the  plot,  representing  the  bit  error  rate  (BER) 
for  each  SNR  level.  The  calculated  error  rate  was  the  total 
number  of  errors  divided  by  the  total  number  of  bits  sent. 
The  results  of  the  simulation  matched  the  general  trend  of 
the  reference  curve,  indicating  an  accurate  model. 


Figure  4 .  Reference  Receiver  Performance  in  AWGN  Over  a 
Range  of  SNR  Values  and  Compared  to  BPSK 
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2. 


Transmitter  Simulation  Design 


The  reference  transmitter  simulation  generated  a  single 
user' s  Bernoulli  binary  random  data  stream  at  the  bit  rate 
and  a  15-stage  shift  register  maximal  length  PN  sequence  at 
the  chip  rate.  The  simulation  mapped  unipolar  binary  out¬ 
puts  to  antipodal  binary  before  mixing.  The  resulting 
chipped  signal  mixed  with  the  carrier  wave,  a  cosine  func¬ 
tion  block  with  a  carrier  frequency  of  1.2288  MHz.  The  re¬ 
sulting  signal  then  passed  through  a  fixed  gain  before  en¬ 
tering  the  AWGN  channel  block.  The  simulation's  parameters 
are  listed  in  Table  1  and  are  based  on  IS-95  parameters  [5] . 


Table  1.  List  of  Reference  Transmitter  Simulation 

Parameters 


Parameter 

Set  Value 

Chip  Frequency  Rc 

1.2288  MHZ 

Bit  Frequency  Rh 

19.2  KHz 

Carrier  Frequency  fc 

4.915200  MHz 

Sample  Time  Ts 

1.27E-08  seconds 

User  Gain  A 

1  (no  units) 

Simulink' s  AWGN  block  provided  the  AWGN  channel  simula¬ 
tion  using  the  parameters  shown  in  Table  2,  which  the  MATLAB 
file  set  prior  to  starting  Simulink. 
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Table  2.  List  of  AWGN  Channel  Simulation  Parameters 


Parameter 

Set  Value 

Initial  Seed 

45  (default  setting) 

Mode 

Signal  to  noise  ratio 

Eb/N0  (dB) 

0  to  9  dB  (set  by  MATLAB 

file) 

Number  of  bits  per  symbol 

1 

Symbol  Period 

A 

The  original  data  goes  to  the  workspace  for  use  in  the 
receiver  and  error  calculations.  The  combined  signal  then 
mixes  with  a  cosine  wave  function  with  a  carrier  frequency 
that  is  four  times  the  chip  rate  and  sampling  frequency  that 
is  64  times  the  chip  frequency.  The  user  gain  is  unitary 
and  has  no  effect  on  the  model  with  the  parameters  as  set. 
Integer  relations  between  the  chip  and  bit  frequencies  were 
arbitrary,  but  we  require  fc>Rc  and  Rs»fc. 

3 .  Reference  Receiver  Model 

The  receiver  shown  in  Figure  5  simulates  the  reference 
receiver  in  Figure  2.  It  consists  of  a  downconverter  and 
filter  subsystem,  a  baseband  processor,  and  an  integration 
block  that  compares  the  demodulated  signal  to  a  threshold  to 
determine  individual  bits.  The  error  rate  calculation  block 
compares  the  threshold  output  against  the  original  data  gen¬ 
erated  by  the  Bernoulli  binary  sequence  block  in  the  trans¬ 
mitter.  However,  this  visual  output  is  for  the  user  and  a 
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separate  sink  saves  each  bit  to  the  workspace  in  a  single 
vector  for  later  error  calculations  and  plotting. 


Original  Data  In 


Convert 


J  Original  Data 


Data  Type  Conversion 


Error  Rate 
Calculation 

T*  Error  Rate 
Rx  Calculation 


Downconvert 
and  Filter 


TB 

Scope 


Reference_RCVR_BER 

r(t)  S*gout — ^Signal  In  DecsVarj — ►  < 

=  Integrate 

nThres 

told  '  - 1 

To  Workspace 

Baseband 

ComP 

are 

Data Recovered 

T o  Workspace 


Figure  5.  SIMULINK  Reference  Receiver  Model 


The  downcoverter  and  filter  subsystem  mixes  the  signal 
with  a  cosine  wave  with  the  same  parameters  as  in  the  trans¬ 
mitter.  Next,  the  signal  passes  through  a  discrete-time 
low-pass  filter  (LPF)  using  a  unitary  pulse  that  is  64  sam¬ 
ples  in  length  to  shape  the  output.  This  simulates  the  ef¬ 
fect  of  the  matched  filter  and  eliminates  the  double  fre¬ 
quency.  Finally,  the  receiver  downsamples  the  output  with 
downsampling  factor  of  64  before  proceeding  to  the  baseband 
processor,  which  simulates  the  reference  receiver's  sampler. 

The  baseband  processor  models  the  de-spreading,  matched 
filter,  and  sampler  from  Figure  2,  and  mixes  the  received 
signal  with  a  PN  sequence  that  matches  the  transmitter.  The 
resulting  signal  passes  through  a  matched  filter  of  64  ones 
and  downsamples  64  times.  The  processed  signal  output  then 
goes  to  the  decision  block  to  determine  bits.  The  integra¬ 
tion  block  performs  the  threshold  comparison  and  decision 
functions  producing  a  unipolar  binary  signal  output,  and 
sends  the  results  to  the  workspace  for  comparison.  The  re¬ 
ceiver  simulation  matched  the  theoretical  performance,  so  it 
provides  a  point  of  reference  for  performance  comparison 
with  the  cyclic-PN  code  receiver  in  the  next  chapter. 
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IV.  CYCLIC  PSEUDO-NOISE  CODE  RECEIVER 


The  cyclic  PN  code  receiver  solves  the  scalability 
problem  while  maintaining  the  BPSK  error  performance  of  the 
reference  receiver.  Its  use  of  DFTs  implemented  through  the 
FFT  algorithm  eliminates  the  need  for  multiple  matched  fil¬ 
ters  to  extract  each  transmitter's  data.  Rather,  a  single 
receiver  block  performs  all  the  same  functions  for  all  the 
transmitters  except  for  bit  decisions. 

A.  CYCLIC  RECEIVER  DESIGN 

The  cyclic  PN-code  design  uses  two  spreading  codes, 
with  each  transmitter  using  a  cyclically  rotated  version  of 
the  same  short  code  mixed  with  a  longer  common  PN  code.  The 
shifted  code  repeats  for  each  bit  period  and  allows  the  re¬ 
ceiver  to  identify  the  transmitter  and  demodulate  data. 
When  mixed  with  the  long  PN  code,  whose  duration  is  much 
greater  than  a  single  bit,  the  resulting  pseudo-noise  se¬ 
quence  is  unique  and  improves  the  LPE  properties  of  the  sys¬ 
tem.  The  DFTs  allow  a  single  receiver  to  analyze  all  trans¬ 
mitted  signals  simultaneously  and  determines  the  number  of 
shifts  of  the  short  PN  code  in  each  signal  as  well  as  the 
data  content.  From  this,  the  receiver  determines  transmit¬ 
ter  identity  and  makes  bit  decisions.  The  theoretical  bit 
error  rate  performance  of  this  receiver  is  the  same  as  the 

f 

reference  receiver,  Ph-Q 

V 
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1 .  Cyclic  PN-Code  Transmitter 


The  cyclically  shifted  PN  code  transmitter  mixes  the 
long  PN  code,  b(t)  ,  with  a  cyclically  shifted  shorter  code, 

c(hf(/?))  ,  where  t  is  time,  y  is  the  length  of  the  shift, 

and  p  is  the  transmitter  number.  The  transmitter  number 
and  quantity  of  shift  are  related  by  an  integer  factor  k  , 
resulting  in  a  relationship  of  y  =  pk  .  For  this  analysis, 
k= 1,  so  the  shift  number  will  be  the  same  as  the  transmit¬ 
ter's  number,  y  =  p.  The  shifts  in  c(V, /(/>))  are  not  calcu¬ 
lated  by  each  transmitter  but  are  preprogrammed  prior  to 
system  emplacement.  The  block  diagram  for  the  transmitter 
is  shown  in  Figure  6. 


The  cyclically  shifted  PN  code  repeats  within  each  bit, 
allowing  for  demodulation  and  identification.  It  can  be 

JV— 1 

represented  as  c(t, y{p\\  =  'y'c\i  —  y{p)\  pTc{ where  p  is  the 

i= 0 

transmitter  number,  t  is  time,  y{p)  is  the  number  of  shifts 
of  the  reference  PN  code  for  transmitter  p  ,  and  PTc(t)  =  1  if 

The  sequence  is  ob- 


0  <t<Tc  and  pTc(t)  =  0  otherwise. 


tained  by  cyclically  shifting  the  N  chip  sequence,  c\i\  ,  y 
shifts  to  the  right.  Each  c[/]e{-l,l}  and  p  does  not  have  to 
equal  y  . 

The  longer  PN  code,  b[t)  ,  provides  security  by  not  re¬ 
peating  for  the  duration  of  the  message.  The  two  PN  codes 
are  generated  and  mixed  at  the  chip  rate  Rc  .  Using  two 
codes  creates  a  unique  spreading  sequence  for  each  transmit¬ 
ter  . 

The  resulting  unique  PN  code  then  mixes  with  the  an¬ 
tipodal  data  signal  mp  before  upconversion  to  a  BPSK  trans¬ 
mission  at  a  predetermined  carrier  frequency  fc  .  The  an- 

oo 

tipodal  data  signal  can  be  represented  as  mp  {t)  =  YJdlPrb  (t~iTb)  , 

-00 

where  df  is  the  ith  bit  and  Pn{t-iTb )  is  the  pulse  function 

for  the  i"‘  bit  time.  Finally,  the  transmitter  amplifies  the 
signal  by  a  pre-determined  factor  A  before  transmission, 
yielding  the  transmitted  signal 

5(7)  =  Ab(t)c  (t,y)  cos  [Ik fct)  .  (4.1) 

2 .  Cyclic  PN  Code  Receiver 

The  signal  in  equation  (4.1)  travels  through  an  addi¬ 
tive  white  Gaussian  noise  (AWGN)  channel  and  arrives  at  the 
receiver  depicted  in  Figure  7  with  the  form  shown  in  equa¬ 
tion  (4.2)  . 

r(t)  =  AYJdlc(t,y)b(t)  cos  (in  ft)  +  «(0  •  (4.2) 

The  signal  enters  the  downcoverter  block  and  changes  to 
the  form  shown  in  equation  (4.3) . 
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a(0  =  djPr  {t-iTbYp  (t,y)b(t)cos2  (27rfct)  +  n(t)cos(27rfct)  (4.3) 


/(«) 


Based  on  the  trigonometric  identity,  cos^  (x)  =  — (l-cos(2x)) 


[11],  this  changes  to  the  form  shown  in  equation  (4.4),  with 
the  signal  received  now  in  three  parts  -  the  desired  signal 
term,  the  double  frequency  term,  and  the  noise  term. 


A  *  A  ^ 

,A_i;L,diPT(t-iTb)F(t>r)b(t)  +^2ldiPr(t-iTby(t’r)b(t)cos(4xfct) 

dll)  A  _oo  A  _oq  \  hi  •  ni  ) 

+ n  (t)  cos  (2  7i  fct) 

Next,  examining  the  desired  signal  term  allows  deriva¬ 
tion  of  the  mean  of  the  decision  statistic.  Noise  will  be 
examined  separately  in  the  next  subsection.  Mixing  with  the 
long  PN  code  effectively  cancels  the  factor  b(t )  from  the 


first  two  terms  on  the  right  hand  side  of  equation  (4.4)  . 
The  LPF  then  removes  the  double  frequency  term,  giving  the 
output 

A  +°°  y-i 

/W=^Zrf.Zc«A[('-P+1K)/7:]  n.5) 

A  /=— oo  7=0 

where  A(x)  =  l  — |x|  if  |x|  <  1  and  A(x)  =  0  otherwise. 
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The  receiver  then  samples  the  signal  f(t)  at  the  chip 

rate  prior  to  it  entering  the  DFT  in  order  to  analyze  a  sin¬ 
gle  bit. 

A 

s[n\  =  —dic[n~r]N  (4-6) 


Analysis  of  the  signal-only  term  leads  to  the  mean  of 
the  decision  statistic.  Starting  from  the  output  of  the  in¬ 
verse  DFT  (IDFT),  the  output  of  the  receiver  is  shown  in  eq¬ 
uation  (4.7).  C(&,0)  represents  the  complex  conjugate  of 

the  DFT  of  the  short  PN  code  with  no  shifts.  Pre¬ 
calculating  this  result  and  storing  it  in  memory  rather  than 
performing  a  real-time  DFT  and  complex  multiplication  helps 
improve  scalability. 


l[n\  =  IDFT  {C(jfc,0)  )DFT[g[n]) 


(4.7) 


Multiplying  the  two  DFT  outputs  is  equivalent  to  circu¬ 
lar  convolution  in  the  time  domain  [12]  .  The  operator  0 
indicates  a  circular  convolution. 


H(k)X(k)  =  DFT{h[n\®x[n]}  (4.8) 

The  complex  conjugate,  indicated  by  * ,  in  the  fre¬ 
quency  domain  is  equivalent  to  the  time-reversed  signal  in 
the  time  domain  [12] .  In  this  instance,  the  complex  conju¬ 
gate  for  the  discrete  time  domain  signal  in  reverse  is  real, 
so  the  complex  notation  drops  out. 

x\n]<^X\k] 

(4.9) 

x[-n\^r[k\ 

A  circular  convolution  is  the  same  as  an  ordinary  dis¬ 
crete  time  convolution,  except  that  the  time-shifted  signal 
is  circularly  shifted,  vice  an  ordinary  shift  [12] . 
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iV  —1 

h  [ft]  0  x  [ft]  =  ^  h  [/ft]  x  [/?  -  /// 


(4.10) 


Applying  equations  (4.8),  (4.9),  and  (4.10)  to  equation 

r  AdiNSn 

(4.7)  yields  an  IDFT  output  of  /[//]  = - - — —  as  shown  in 

equation  (4.11)  .  This  holds  for  0  <n  <  N  —  1  ,  where  Sn  is  the 
Kronecker  delta  and  it  is  assumed  that  c[//]  has  low  correla¬ 
tion  with  c[n-i\N  for  all  0<i<N-l.  The  square  of  a  shifted 
PN  codes  is  one  and  the  sum  of  that  result  over  /V -terms  is 


[n]  =  IDFT^\c(k,^]DFT(g[n]) 

=  &  [«]  ®  c  [—ft] 

=  yc[”-40cH 

Ad  £4 

=  y1Ec[w-^Lc[/ft-ft]JV 

^  m= 0 

AdiNSn/ 


(4.11) 


When  n^y  ,  using  a  circular  convolution  definition 
with  a  Kronecker  delta  gives  a  result  of  0.  Choos- 

JV-l 

ing  Tc\m-y\,c\m-n\^  NSyn  allows  the  result  shown  in  equa- 

m= 0 

tion  (4.11) .  In  the  case  of  multiple  transmitters'  signals 
arriving  simultaneously,  /[//]  indicates  whether  transmitter 

ft  is  transmitting  (if  |/[/?.]|  is  large)  and  the  value  of  the 
sent  data  bit  (sgn(/[ft]j).  Thus  / [ft]  for  0<n<N-l  simultane¬ 
ously  indicates  which  of  the  N  transmitters  are  active  and 
the  bits  they  have  transmitted. 
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3. 


Noise  Calculations 


Noise  at  r(t )  is  additive, 
ju  =  0  ,  <j2  =  +oo ,  and  a  PSD  of 

Nn 

changes  the  PSD  to  due  to 

squared  cosine  function,  which  is 


white 

2 

the 

1/ 

72  • 


and  Gaussian  and  has 

The  downconversion 


average  power  of  the 


After  the  LPF,  the  variance  of  the  noise  is  no  longer 
infinite,  but  limited  by  the  bandwidth  of  the  filter. 


(4 . 12) 


N0RhN 


where  //(/)  is  the  frequency  response  of  the  LPF  and  sam¬ 
pling  does  not  affect  white  Gaussian  noise.  Therefore,  the 
noise  at  g[»]  ,  the  input  to  the  DFT,  is  still  a  Gaussian 


random  variable  with  zero  mean 


and 


2  N0RbN 
(j  =  — — - —  ,  or 


( 


N 


n  W 
’  4 


The  samples,  g[»]  ,  are  independent  because  the 


input  to  the  LPF  is  white  noise,  the  LPF' s  impulse  response 
has  a  duration  equal  to  the  chip  duration,  and  the  LPF' s 
sample  rate  is  equal  to  the  chip  rate.  Therefore: 


E{g[m\s[p]}  =  N{)RhNA5"lp 


(4 . 13) 
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The  output  of  the  mixer  in  the  frequency  domain  is: 


L(k)  =  G(k)C(k  |0)*.  (4.14) 

The  expected  output  of  the  receiver  is  then  the  IDFT  of 
L{k ) .  The  product  of  the  two  functions  in  the  frequency  do¬ 
main  is  equivalent  to  a  circular  convolution  in  the  time  do¬ 
main  [  12  ]  . 

l[n]  =  IDFT\G(k)C{k\Q)*} 

=  g[«]®c[-n]  (4.15) 

N- 1 

=  J]s[m\c[m-n\N 

m= 0 

The  expected  value  or  mean  of  the  noise  component  of 
this  function  is: 

E  M}  =  E  (X  S  M  c [m  -  n]N  | 

Lm=0  J 

=  YjE{s[m\}c[m-n\N  (4.16) 

m= 0 

=  0 

The  variance  of  this  function  is  the  variance  of  the 
receiver  output  and  therefore  the  variance  of  the  signal 
(which  allows  probability  of  error  calculations) , 

{JV-l  JV-l 

Z  s  [™)c  [m -n)NYJg{p)c{p-  n)N 

m= 0  p= 0 

N-l  N- 1 

=  THE{g(m)g{p)}c(m-n)Nc{P-n)N 

m= 0 p= 0 

N-l  ,  .  2 

=  Tj E\ g{m)  }  [c(m  ~n)N]  (4.17) 

m= 0  V  ^ 

m= 0  4 

KRhN2 

4 
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4. 


Theoretical  Error  Performance 


At  the  output  of  the  IDFT  block,  the  signal's  decision 

This  is  used  to  find  the 


.  .  .  ,r  ,  i  Ad.N  N0RhN 

statistic  is  l[n]~N'  h 


V 


J 


probability  of  bit  error,  Q 


2E, 


K'l  N°  J 


,  the  same  performance  as 


a  BPSK  system.  Using  the  mean  and  variance  calculated  above 
in  equations  (4.11)  and  (4.17),  and  the  equivalence 
Eh  =  0.5A2Th,  the  bit  error  rate  is 

Pb  =  Pr(/(«)  >  0  |  dj  =  -l) 

AN  AN 


=  Pr 


/ («)+- 


N0RbN 


> 


=  Q 


2  E, 


4 

-1 


V 


N0RbNz 


(4.18) 


\\  No  J 


B. 


MODELING  THE  CYCLIC  PN-CODE  RECEIVER 


The  proposed  receiver  was  simulated  in  MATLAB  and 
tested  for  100,000  bits  at  each  of  a  range  of  signal  to 

from  0  dB  to  9  dB .  A  single  transmit¬ 
ter  was  simulated  with  a  one  chip  cyclic  shift  (/= 1)  in  the 
base  short  PN-code.  The  simulation  used  an  AWGN  channel  and 
BPSK  modulation  for  the  carrier.  Perfect  phase  synchroniza¬ 
tion  was  assumed.  The  MATLAB  code  for  the  simulation  is  in 
Appendix  B. 


noise  ratios 
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1. 


Simulation  Results 


The  cyclic  PN  code  receiver  simulation  results  are 
shown  in  Figure  8 .  The  black  line  represents  the  BPSK  theo¬ 
retical  bit  error  rate,  and  is  the  same  curve  from  Figure  3. 
The  simulation  bit  error  rate  results  are  the  blue  squares 
on  the  plot,  representing  the  BER  for  each  SNR  level.  As 
with  the  reference  receiver,  the  calculated  error  rate  was 
the  total  number  of  errors  divided  by  the  total  number  of 
bits  sent.  The  results  of  the  simulation  matched  the  gen¬ 
eral  trend  of  the  analytically  derived  curve,  indicating 
that  model's  accuracy. 

Figure  9  is  an  example  of  the  IDFT  output  for  a  single 
bit  from  a  single  transmitter  with  a  single  shift  of  the 
reference  short  PN  code,  i.e.,  y  =  1.  The  signal  from  the 
transmitter  of  interest  stands  out  clearly  from  the  64  pos¬ 
sible  signals  received  during  this  bit  time. 


Figure  8.  Cyclic  PN  Code  Receiver  Performance  in  AWGN  over  a 

Range  of  SNR 
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IDFT  output  for  bit  #1 ,  SNR=9  dB 


Figure  9.  IDFT  output  for  a  single  bit  duration,  at 


=  9  dB  with  y  =  1 


Figure  8  shows  good  concurrence  between  simulation  and 
theory.  Figure  9  shows  that  lini  at  the  receiver  indicates 

reception  of  a  +1  from  transmitter  number  one  and  only  noise 
from  the  other  transmitters.  A  mechanism  is  still  needed  to 
distinguish  an  active  transmitter's  signal  from  noise  only. 
This  is  not  addressed  in  this  thesis  except  to  note  that 
conventional  methods  such  as  using  start  and  stop  bit  pat¬ 
terns  or  employing  two  thresholds  to  make  "+1",  "not  trans¬ 
mitted"  and  "-1"  can  be  used  to  address  this  issue  [13] . 

2.  Transmitter  Simulation  Design 

The  simulation  of  the  cyclic  PN  code  transmitter  used  a 
fixed  i?6=l  and  then  generated  other  variables  from  its  mul- 
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tiples.  Unlike  the  reference  receiver,  these  differences 
were  relative  and  not  on  an  absolute  scale.  The  simulated 
transmitter  mixed  the  fixed  data  sequence  with  two  separate 
PN  codes,  one  of  which  was  the  cyclic  rotation  of  the  short 
PN  code.  The  resulting  signal  was  sampled  at  the  chip  rate 
Rc  =  64Rh  and  then  upsampled  with  upsampling  factor  of  256  be¬ 
fore  mixing  with  the  cosine  vector  at  frequency, 
fc=16Rc=\024Rb,  and  a  fixed  amplitude.  The  transmitter  pa¬ 
rameters  are  in  0 . 


Table  3.  Parameters  for  Cyclic  PN  Code  Transmitter 

Simulation 


Parameter 

Value 

Rh  (bit  rate) 

1 

Rc  (chip  rate) 

64  Rh 

fc  (carrier  frequency) 

\6Rc=\024Rb 

i?s  (modulated  signal  sample  rate) 

16/c 

Ts  (sample  time) 

A 

Data  bits 

50,000 

The  cyclic  PN  code  simulation  did  not  use  the  AWGN 
block  in  SIMULINK.  Rather  it  generated  noise  as  a  vector  of 
the  same  length  as  a  single  upsampled  bit  based  on  the  SNR 
for  a  particular  loop  and  the  noise  PSD  using 
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Noise= 


R  N 

s  o 

2 


randn(l,Rs)  . 


The  noise  PSD  was  then  based  on  the 


i  j  , 

energy  per  bit  by  Na  = — ,  with  Eh  fixed  at  0.5.  The 

10'  /l0 

MATLAB  function  "randn"  generated  a  vector  of  random  numbers 
with  length  equal  to  the  sampling  rate  with  a  coefficient 
based  on  the  noise  PSD,  allowing  variable  noise  based  on  SNR 
and  vector  addition.  The  effect  on  the  transmitted  signal 
was  equal  to  the  effect  of  AWGN. 


3 .  Cyclic  PN  Code  Receiver  Model 


The  cyclic  PN  code  simulation  downconverted  the  re¬ 
ceived  signal  with  noise  by  multiplying  the  data  vector  by 
the  same  cosine  vector  used  in  the  transmitter.  To  de¬ 
spread  the  signal  with  the  long  PN  code,  the  simulation  up- 
sampled  the  PN  code  before  mixing.  A  square  pulse  whose 
length  was  the  chip  duration  performed  the  low  pass  filter 
function  through  convolution  with  the  de-mixed  signal.  The 
simulation  sampled  the  LPF  output  at  the  chipping  rate,  by 
downsampling  using  the  ratio  between  sampling  and  chipping 
rates . 

MATLAB' s  Fastest  Fourier  Transform  in  the  West  (FFTW) 
algorithm,  the  "fft"  function,  acted  as  the  DFT  block  pro¬ 
posed  in  the  theoretical  model.  The  simulation  directly 
calculated  the  DFTs  of  the  signal  received  and  reference 
short  PN  code,  the  complex  conjugate  of  the  short  PN  code, 
and  mixed  the  DFT  outputs.  The  MATLAB  function  "if ft"  acted 
as  the  IDFT  block,  and  a  loop  made  bit  decisions  based  on 
the  IDFT  output. 
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The  cyclic-PN  code  receiver  simulation  performed  well 
and  matched  the  theoretical  BPSK  error  rate.  In  addition  to 
the  cyclic-PN  code  receiver,  a  coding  mask  receiver  was  pro¬ 
posed  as  a  possible  alternative.  The  next  chapter  analyzes 
the  coding  mask  receiver  and  compares  its  theoretical  error 
performance  to  the  reference  receiver. 
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V.  CODING  MASK  RECEIVER 


As  an  additional  point  of  comparison,  one  possible  al¬ 
ternative  to  the  cyclic  PN  code  receiver  is  a  coding  mask 
receiver.  This  design  relies  on  similarities  in  the  trans¬ 
mitters'’  PN  codes  to  perform  a  two-stage  detection  and  de¬ 
modulation.  This  receiver  is  not  a  viable  solution  as  a 
system  if  many  transmitters  are  expected  to  be  simultane¬ 
ously  active  as  it  does  not  improve  scalability  in  this  case 
and  it  has  a  greater  error  rate  than  the  reference  receiver. 
However,  it  requires  fewer  processes  for  single  transmitters 
and  could  be  a  viable  alternative  design. 

A.  RECEIVER  DESIGN 

The  coding  mask  receiver  subdivides  the  full  body  of 
transmitters  into  groups,  which  will  be  called  subsets  for 
the  remainder  of  this  thesis.  Each  subset  has  related  PN 
codes  with  some  common  chips  to  speed-up  processing  and  this 
is  the  "mask."  The  receiver's  signal  detection  circuit  mix¬ 
es  the  combined  signal  received  with  the  base  PN  code  of  its 
subset  and  then  its  mask  to  determine  signal  presence.  When 
a  signal  from  a  particular  subset  is  detected,  all  receivers 
within  that  subset  process  in  parallel  to  extract  the  signal 
received.  This  design  has  non-linear  growth  and  while  it  is 
of  similar  complexity,  efficiency  is  gained  through  savings 
in  total  number  of  processes.  This  savings  would  be  most 
beneficial  in  a  software  defined  radio  implementation  since 
the  savings  in  processing  is  realized  by  simpler  and  faster 
algorithms  and  reduced  requirements  on  the  microprocessor. 
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digital  signal  processor,  or  field  programmable  gate  array 
(FPGA)  responsible  for  the  signal  processing  [14] . 

The  coding  mask  receiver  system  would  use  the  same 
transmitter  as  the  reference  receiver  as  shown  in  Figure  1 
and  generate  a  similar  transmitted  signal  in  the  form 
s(t)  =  Amp  (t)c  (t)cos(2?rfct)  .  Here,  mp(t)  is  the  antipodal  data 
with  the  subscript  p  indicating  the  transmitter  number  and 
having  the  same  definition  as  equation  (3.1) .  The  PN  code, 
c  (t) ,  is  unique  and  antipodal  for  each  transmitter  and  the 
subscript  again  indicates  the  transmitter  number.  The  local 
oscillator,  co%{lKfct)  ,  is  also  BPSK,  and  the  coding  mask 
transmitter  also  amplifies  the  transmitted  signal  by  a  pre¬ 
set  gain  A  . 

The  signal  then  passes  through  an  additive  white  Gaus¬ 
sian  noise  (AWGN)  channel,  which  adds  the  noise  term  n(t )  to 
the  combined  signal.  The  signal  received  is,  therefore, 
r(t}  =  Amp  (V)c  {t}cos(27ifct}  +  n{t}  .  This  signal  then  enters  the 
coding  mask  receiver  shown  in  Figure  10. 


38 


Figure  10.  Code  Masking  Receiver  Diagram 


In  Figure  10,  a  receiver  for  a  single  subset  is  shown. 
The  PN  codes  of  this  subset  have  some  common  points  allowing 
the  mask  to  be  pre-calculated  and  preprogrammed  into  the  de¬ 
tection  branch.  The  first  PN  code  in  this  subset  is  c0(V)  , 

where  the  subscript  0  indicates  the  ordinal  numbering  of 
the  transmitters  in  the  subset.  The  first  stage  accom¬ 
plishes  signal  detection  by  mixing  the  downconverted  signal 
with  the  first  PN  code  in  the  subset  c0(V),  the  carrier 

cos(2^/c?) ,  and  the  mask  k[t )  before  entering  the  matched  fil¬ 
ter.  The  mask  represents  the  common  points  of  all  PN  codes 
in  a  subset,  with  1  used  to  indicate  they  are  the  same  and 
0  used  if  even  a  signal  code  is  different  at  a  particular 
chip.  The  results  is  £(V)e{l,0{  and  not  e{l,-l}.  The  signal 

received  at  this  point  is  similar  as  in  the  reference  re¬ 
ceiver  with  three  components,  a  signal  of  interest,  a  double 
frequency,  and  noise  [7] . 
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A  A 

r'{t)  =  —m{t)k(t)-  —  m(t)k(t)cos(Anfct)  +  n{t)cQ(t)k{t)cos(2nfct)  (5.1) 
The  matched  filter's  impulse  response  h(t)  is  a  pulse 


function  pn(t-iTh^ ,  where  the  value  of  the  function  is  1 

during  the  bit  interval  and  0  otherwise.  The  resulting  out¬ 
put  is  the  signal  component  of  interest  output  [7] . 


y(<) 


A 

=  —  m 
2 


(t)k(t)*h(t) 
m(a)k(a)h(t  -a)  da 


(5.2) 


However,  the  integrand  is  only  non-zero  when  t  —  Tb<a<t, 
changing  the  integration  interval.  A  substitution  can  be 
made  for  the  modulated  signal,  whereby  the  antipodal  data, 
m(tf,  can  be  represented  as  a  single  bit  at  a  discrete  time 
instance  rather  than  as  the  sum  of  bits  and  the  pulse  func¬ 
tion,  as  was  shown  earlier:  m(t)  =  J'idipT(a—iTh)  =  df_ ,  .  Addi¬ 
tionally,  the  impulse  response  of  the  matched  filter  taken 
at  a  discrete  time  instance  where  t  =  tTh  results  in  the  sub¬ 
stitution  h(lTb-a)  =  pT{lTb-a)  .  Using  these  two  substitutions 

changes,  the  matched  filter  output  in  the  detection  branch 
to  a  different  form  of  the  mean  of  the  signal  received  at 
that  point. 

t 

J  m(a)k(a)h(t -a)da 

t-rb 

A 

y(lTb)  =  —  j  d^kioc)  pT(dTh-a)da  (5.3) 

2  U-l  )T„ 

=  y  df-\PTh 
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This  result  is  identical  to  the  reference  receiver' s 
signal  of  interest  with  the  exception  of  /? ,  the  fraction  of 

Th  for  which  the  mask  is  one,  i.e.,  (3  =  —  J  k{t)dt  .  Since  the 

o 

matched  filter  is  a  low  pass  filter,  the  double  frequency 
term  is  filtered  out. 


N 

The  PSD  of  the  noise  at  the  receive  antenna  is  S  =  — -. 

2 

As  in  the  reference  receiver,  when  this  noise  is  mixed  with 

Nn 


the  carrier,  the  resulting  PSD  is  S  =  — 

4 


However,  the  mask 

has  two  possible  states,  k[t)  =  0  or  *(0=i  •  In  the  former 

case,  the  resulting  mean  { ju  )  ,  variance  ( a2  )  ,  and  PSD  ( S  ) 
of  the  noise  would  be  zero.  In  the  latter  case,  //  =  0, 

cr2=+co,  and  S  =  ^~,  which  is  white  Gaussian  noise.  However, 

4 

to  account  for  the  fraction  of  the  mask  that  equals  1  in  a 
given  period  of  Tb  ,  the  term  /?  can  likewise  be  used  as  in 
the  signal  of  interest,  giving  a  resulting  input  PSD  of 
N  B 

S  = — - — .  The  noise  PSD  at  the  output  of  the  matched  filter 

m  4 


is  therefore  Sout(f)  =  \H  (f) 
power  at  the 


■V 


Its  integral  is  the  noise 


matched  filter  output 

H(f)\  df  .  Applying  Parseval's  theo¬ 


rem  [12],  changes  this  to  o'  =  J  h 1  ( t)dt  .  And  this  in  turn 

4  A 


Tb 


solves  to  the  noise  power  or  variance 
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_  n0/3t„ 

4 


PERFORMANCE  ANALYSIS 


Using  the  sampled  matched  filter  output  for  the  signal 
of  interest  and  noise,  the  decision  statistic  for  signal  de¬ 


tection  is  N 


Adt_xprb  N0(3Th 


The  probability  of  bit  error 


is  the  same  as  the  probability  of  deciding  a  "1"  was  sent 
when  a  "-1"  was  sent.  As  shown  in  equation  (5.5),  the  prob¬ 
ability  of  error  for  a  "-1"  sent  is  equivalent  to  the  prob¬ 
ability  that  the  mean  of  the  signal  plus  the  noise  N  is 
greater  than  zero. 


Pr  (4  =  1 1  d,_t  =  -1)  =  Prf  -^  +  N  >  0 

V  2 


2  Ebp 


(5.5) 


In  terms  of  signal  power  required,  the  best  possible 
condition  would  be  /?  =  1 ,  when  all  PN  codes  are  the  same. 
However,  this  condition  would  make  multiple  access  impossi¬ 
ble— the  receiver  would  be  unable  to  distinguish  between  the 
signals . 
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BER 


Figure  11.  BPSK  Error  Performance  for  Fixed  SNR 

and  variable  /? 
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Effects  of  Beta  on  BER 


Figure  12.  Comparison  of  Theoretical  BPSK  and 
Three  Values  of  /? 


Fewer  chips  in  common  only  increases  the  bit  energy  re¬ 
quired  to  achieve  a  similar  bit  error  rate  as  the  reference 
receiver.  The  resulting  system  probability  of  bit  error  is 
worse  than  reference  receiver  for  the  same  signal  to  noise 
ratio.  The  first  Q  function  is  the  error  of  the  detection 
branch  where  /?  represents  the  effect  of  the  mask.  The  sec¬ 
ond  Q  function  is  the  demodulation  error  of  the  receiver 
branches  in  the  subset.  The  result  is  a  performance  worse 


than  2 Q 


2  E„ 


N, 


o  J 


p,  M  =  Q 


2  Eb0 


N, 


+  Q 


0  J 


2  Ek 


N, 


(5.6) 


o  J 
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This  receiver  was  not  simulated  due  to  time  constraints 
and  its  significantly  poorer  BER  performance  as  compared  to 
the  receiver  in  Chapter  IV.  Additionally  because  of  its  in¬ 
ferior  theoretical  error  performance,  the  overall  system 
performance  was  not  modeled  and  compared  to  the  reference  or 
cyclic  PN  code  receivers.  Chapter  VI,  therefore,  does  not 
include  the  mask  receiver' s  performance  in  its  comparison 
and  analysis. 
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VI.  RESULTS  AND  COMPARISONS 


The  reference  receiver  and  cyclic  PN  code  simulations 
ran  were  of  the  same  duration,  100,000  bits  per  dB  level. 
Their  results  both  closely  matched  the  theoretical  un-coded 
BPSK  performance  in  an  AWGN  channel,  as  expected.  Of  note, 
both  simulations  did  not  address  fading  or  synchronization 
issues . 

A.  BIT  ERROR  RATE  PERFORMANCE 

The  reference  receiver  and  cyclic  PN  code  DFT  receiver 
have  identical  BER  performance,  as  shown  in  the  analysis  in 
Chapters  III  and  IV,  and  confirmed  in  the  simulation  results 
shown  in  Figure  13.  In  both  cases,  a  SNR  of  approximately  7 
dB  was  required  to  achieve  an  error  rate  of  1  in  1000  and 
approximately  8.5  dB  for  1  in  10000. 
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Figure  13.  Comparison  of  BPSK  Theoretical  (curve).  Reference 
Receiver  (circles)  ,  and  Cyclic  PN  Code  Receiver  (squares) 

B.  SCALABILITY  COMPARISON 

Examining  the  reference  receiver  system  with  N  trans¬ 
mitters,  it  initially  processes  the  received  signals  with  a 
common  antenna,  low-noise  amplifier,  local  oscillator,  and 
low  pass  filter.  From  there,  such  a  receiver  requires  a 
unique  matched  filter,  sampler,  and  bit  decision  block  for 
each  transmitter  in  the  system. 

The  cyclic  PN  code  design  in  Figure  7  and  the  reference 
receiver  in  Figure  2  have  similar  components  from  antenna  to 
LPF.  Looking  at  the  number  of  multiplications  required  al¬ 
lows  for  a  scalability  comparison.  For  the  reference  re¬ 
ceiver,  there  are  N  branches  from  the  LPF  forward.  Each 
branch  must  analyze  N  signals,  with  despreading,  filtering. 
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sampling,  and  bit  decisions.  Therefore,  a  good  estimation 
of  the  scalability  is  total  number  of  multiplications  re¬ 
quired  for  despreading  in  all  the  branches.  Sign  multipli¬ 
cations  (recall  c[«]e{-l,l})  are  not  as  complicated  as  other 

multiplications,  so  a  coefficient,  a<  1,  is  used  to  indicate 
the  ratio  of  the  complexity  of  a  sign  multiplication  to  the 
complexity  of  a  complex  number  multiplication.  For  each 
branch  in  the  linear-growth  receiver,  there  are  N  sign  mul¬ 
tiplications  per  bit  duration  because  there  are  N  chips  per 
bit.  Therefore,  the  total  effective  number  of  complex  mul¬ 
tiplies  per  bit  duration  is 

C reference  =N(a  N) 

(0.1) 

=  N2a 

This  value  is  accepted  as  the  measure  of  complexity  of  the 
reference  receiver. 

For  the  cyclic  PN  code  receiver  operating  in  the  same 
scenario,  it  must  complete  N  sign  multiplications  per  bit 
duration  to  despread  the  long  PN  code,  N  complex  multiplies 
in  the  last  mixer,  1  DFT,  and  1  IDFT.  The  FFT  and  IFFT  each 
require  0.5N\og2N  complex  multiplies  [6].  Combining  these 
gives  the  total  number  of  complex  multiplies  per  bit  dura¬ 
tion  for  the  proposed  receiver,  which  is  accepted  as  its 
measure  of  complexity. 

C^«-r„=ciN  +  y2log2N  +  N  +  y2\og2N 
=  N(\og2N  +  l  +  a) 

The  ratio  of  the  complexity  of  the  reference  receiver 
to  the  complexity  of  the  cyclic  PN  code  receiver  is 
Na 

- ,  where  the  value  of  a  depends  on  the  bit  resolu- 

log2  N  + 1  +  a 
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tion  of  the  numbers  multiplied.  The  complexity  of  one  com¬ 
plex  multiply  with  real  and  imaginary  parts  each  represented 
with  r  bits  is  equivalent  to  the  complexity  of  4 r  sign  mul¬ 
tiplications.  If  the  number  of  transmitters  exceeds  ap¬ 
proximately  300,  the  proposed  receiver  is  a  more  efficient 

implementation,  assuming  a  =  ^2  '  w^1-'-c^L  corresponds  to  eight 

bit  resolution.  The  relative  growth  in  the  size  of  the  two 
systems  is  shown  in  Figure  14. 


x  io4  Comparing  Receiver  Scalability 


Figure  14 . 


Scalability  Comparison  for  7^  =  1000 


and 
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VII .  CONCLUSIONS  AND  RECOMMENDATIONS 


As  can  be  seen  in  Figure  8  and  Figure  13,  the  proposed 
design  matches  theoretical  uncoded  BPSK  performance  (which 
also  matches  the  reference  receiver) .  The  IDFT  output,  as 
shown  in  Figure  9  for  a  single  transmitter,  clearly  indi¬ 
cates  the  shift  number  and  bit  value.  The  simulation  used 
an  N  =  64-chip  PN  code,  allowing  up  to  64  transmitters  in 
the  model  system.  For  a  full  system,  64-bit  decisions  would 
be  required  for  each  Th  interval.  N  =  64  was  chosen  for  il¬ 
lustration  purposes  only,  and  much  larger  values  are  practi¬ 
cal.  The  cyclic-PN  code  receiver  maintains  the  BER  while 
improving  scalability,  making  it  a  viable  solution. 

A.  PERFORMANCE  COMPARISON 

1 .  Reference  vs .  Cyclic  PN  Code 

Since  the  receivers  had  similar  error  performance  for 
an  uncoded  signal  in  AWGN,  it  stands  that  the  key  difference 
is  receiver  complexity.  In  this,  the  cyclic  PN  code  re¬ 
ceiver  has  superior  scalability  when  the  number  of  transmit¬ 
ters  is  large,  as  shown  in  equations  (6.1)  and  (6.2),  as 
well  as  in  Figure  14. 

2 .  Reference  vs .  Code  Mask 

The  theoretical  bit  error  rate  performance  of  the  code 
mask  receiver  was  inferior  to  the  reference  receiver.  The 
code  mask  receiver' s  advantage  lies  in  the  fewer  required 
processes  to  receive  a  single  signal,  but  its  higher  error 
rate  offsets  this  advantage. 
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3. 


Code  Mask  vs .  Cyclic  PN  Code 


The  cyclic  PN  code  receiver  has  the  same  bit  error  rate 
performance  as  the  reference  receiver,  but  superior  seal- 
ability  and,  therefore,  is  the  better  choice  when  the  number 
of  transmitters  is  large.  Additionally,  the  DFT-based  fil¬ 
tering  in  the  cyclic  PN  code  receiver  was  more  efficient  in 
terms  of  space,  but  potentially  less  efficient  for  a  single 
signal  as  the  total  number  of  processes  required  is  the  same 
regardless  of  the  number  of  transmitters  in  the  system. 
However,  for  a  large  number  of  transmitters,  the  coding  mask 
receiver  is  inefficient,  as  it  requires  two  stages  of  proc¬ 
esses,  resulting  in  more  overall  complexity  than  either  the 
reference  or  cyclic  PN  code  receivers.  The  coding  mask  re¬ 
ceiver'  s  theoretical  BER,  even  under  the  best  conditions 
would  be  twice  that  of  the  cyclic  PN  code  receiver. 

B .  FUTURE  WORK 

1 .  Improving  the  Cyclic  PN  Code  Model 

The  models  presented  were  simulated  in  an  AWGN  environ¬ 
ment.  Simulating  the  same  models  in  a  fading  and  jamming 
environment  would  better  match  potential  real-world  scenar¬ 
ios.  Large  networks  of  sensors  deployed  in  a  mountainous  or 
urban  environment  would  suffer  fading  or  jamming  effects  in 
addition  to  AWGN.  Validation  of  the  cyclic  PN  code  re¬ 
ceiver'  s  error  performance  in  these  environments  is  impor¬ 
tant  . 

MATLAB' s  "randint"  function  and  Simulink's  equivalent 
block  generated  the  simulations'  random  unipolar  binary  se¬ 
quences.  Future  researchers  could  develop  a  truly  orthogo- 
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nal  cyclically  shifted  PN  code  and  validate  its  performance, 
especially  the  prevention  of  inter-transmitter  interference. 
Additionally,  future  research  should  synthesize  an  effective 
forward  error  correction  code  using  BPSK  or  QPSK  for  the  cy¬ 
clic  PN  code  receiver. 

Quantifying  the  LPE  features  expands  the  validity  of 
the  proposed  communications  system.  This  would  require  a 
larger  model  with  multiple  transmitters  and  a  non- 
cooperative  interceptor.  Additionally,  models  of  the  most 
likely  operating  environments,  which  account  for  potential 
adversaries,  are  necessary  in  order  to  analyze  potential 
system  performance. 

This  thesis'  models  also  need  a  synchronization  analy¬ 
sis.  Both  the  reference  and  cyclic  PN  receiver  in  this  the¬ 
sis  relied  on  assumed  synchronization  to  make  bit  decisions. 
Synchronization  blocks  or  a  non-coherent  form  of  modulation 
might  solve  this  problem. 

2 .  Improving  LPI  and  Data  Rate  Performance 

Strong  encryption  would  improve  the  systems  defenses 
against  exploitation  of  the  intercepted  signal.  However, 
encryption  would  not  improve  overall  LPE  as  adversaries 
could  potentially  exploit  signal  presence  and  geolocation 
through  direction  finding  to  localize  transmitters. 

The  systems  designed  and  tested  in  this  thesis  used 
BPSK  modulation.  Modulation  techniques  like  QPSK,  M-PSK,  M- 
frequency  shift  keying  (MFSK) ,  and  M-quadrature  amplitude 
modulation  (MQAM)  were  not  tested.  Other  modulation  tech¬ 
niques  could  achieve  higher  data  rates  with  trade  offs  in 
transmitter  power,  error  rate,  and  exploitation. 
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APPENDIX  A.  REFERENCE  RECEIVER  SIMULATION  CODE 


All  code  in  this  appendix  was  generated  using  MATLAB 
v7 . 6  (R2008a) .  The  SIMULINK  file  that  supports  this  simula¬ 
tion  is  provided  via  electronic  media. 

%  Name:  LCDR  Frank  Cowan 

%  Title:  Reference  Transmitter  and  Receiver 
%  Simulation  Project:  NPS  EE  Thesis 
%  Date  Created:  22  Jan  2009 
%  MATLAB/Simulink  Version:  7.6.0,  R2008a 

%  This  file  provides  the  variable  inputs  for  a  simulation  of  the 
reference  transmitter  and  receiver  for  my  thesis.  All  data  is  saved 
to  files  and  will  be  plotted  using  a  separate  m-file. 

clc;  %  clears  the  command  window 
clear  all;  %  clears  workspace 

%  Transmitter  parameters 

Chip  Freq  =  1.2288e6;  %  Rc,  set  at  1.2288  MHz 
Bit  Freq  =  Chip  Freq/64;  %  Rb,  set  at  19200  Hz,  1/64  of  Rc 
Carrier  Freq  =  4*Chip  Freq;  %  4915200  Hz,  frequency  of  carrier 
CarrierSampleTime  =  1 / ( 1 6*Carrier_Freq) ;  %16  samples/cycle 

TransmitterGain  =  1;  %  A,  for  Arbitrary  gain 

%AWGN  channel  parameters 

EbNo_array  =  [0,  1,  2,  3,  4,  5,  6,  7,  8,  9];  %  the  EbNo  increments 
InputSignalPower  =  (TransmitterGainA2 ) /2 ;  %  used  in  the  AWGN  channel 
block  for  input  signal  power 

AWGN  bps  =1;  %  used  in  the  AWGN  channel  block  for  bps  setting 
AWGN  symbol  period  =  1/Bit  Freq;  %  sets  the  symbol  period  for  the 
AWGN  in  the  channel 

%  Reference  Receiver  parameters 

IntegrationThreshold  =  0;  %  required  parameter  for  the  reference  re¬ 
ceiver  decision  block 

%  Conduct  simulation 

Run_No  =1;  %  sets  the  Run  No  to  default  setting  of  1 
sim_length  =  3;  %  same  simulation  length  used  for  each  loop  -  based 
on  the  time  interval  used,  this  works  out  to  around  50,000  samples 
per  loop 

for  EbNo  idx  =  1:10 
%  Simulation 

Run  No  %  displays  the  run  number  in  the  work  space  to  allow  user 
to  observe  progress 

EbNo  =  EbNo  array (EbNo  idx);  %  assigns  the  value  of  EbNo  used  in 
the  AWGN  block  in  the  channel 
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calls  the  simula- 


sim (' Reference  RCVR  Model  VI',  sim  length);  % 
tion  and  assigns  the  sim  length 

%  Calculate  the  number  of  errors 
Bits  sent  =  length (Data  Transmitted) -3 ;  %  number  of  bits  sent  by 
transmitter.  -3  is  needed  to  account  for  delays  in  the  system. 

Errors  =  sum (xor (Data  Transmitted,  Data  Recovered) ) -1 ;  %  calcu¬ 
lates  the  number  of  bit  errors  by  XOR'ing  the  data  transmitted  and 
data  received. 

Error  Rate (EbNo  idx) = (Errors/Bits  sent);  %  calculates  the  error 
rate  by  dividing  the  number  of  errors  by  the  total  number  of  bits 
sent.  This  is  the  data  that  will  be  plotted  after  the  loop. 

%  save  the  number  of  bits  created  for  the  simulation  run 
Simulation  Length (EbNo  idx)  =  Bits  sent;  %  saves  the  value  of 
Bits-sent  to  an  array  for  later  analysis 

%  increment  counter 

EbNo  idx  =  EbNo  idx  +1;  %  increments  the  EbNo  index  and  walks 
through  the  EbNo  array 

Run  No  =  Run  No  +1;  %  increments  the  counter  for  the  loop 

end 

save ('Error  Rate',  'Error  Rate')  %  Saves  the  Error  Rate  variable  in¬ 
to  a  separate  file 

save ('EbNo  array',  'EbNo  array')  %  Saves  the  EbNo  array  variable  in¬ 
to  a  separate  file 

%  End  of  Reference  RCVR 
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APPENDIX  B.  CYCLIC  PN-CODE  RECEIVER  SIMULATION 

CODE 


All  code  in  this  appendix  was  generated  using  MATLAB 
v7 . 6  (R2  00  8a)  . 

A.  CYCLIC  PN  CODE  RECEIVER  SIMULATION  MATLAB  CODE 

%  Name:  LCDR  Frank  Cowan 
%  Title:  Cyclic  PN  receiver  Code 
%  File  name:  PN_Code_RCVR_VF .m 
%  Simulation  Project:  NPS  EE  Thesis 
%  Date  Created:  06  May  09 
%  MATLAB/Simulink  Version:  7.6.0,  R2008a 

%  Linked  functions:  upsample2.m 

%  This  files  simulates  a  single  transmitter  system  for  the  cyclic  PN 
code 

%  transmitter  and  receiver  system  detailed  in  my  thesis. 

clc;  %  clears  the  command  window 
clear  all;  %  clears  workspace 

%  create  data  array 

Rb  =  1;  %  bit  rate 

Rc  =  64*Rb;  %  chip  rate 

alpha  =  16;  %  ratio  of  fc  to  Rc 

fc  =  alpha*Rc;  %  carrier  frequency  =  alpha  factor  times  Rc 
beta  =  16;  %  ratio  of  Rs  to  fc 
Rs  =  beta*fc;  %  sample  rate 

Ts  =  1/Rs;  %  sampling  time  -  inverse  of  sampling  rate 
N  =  49999;  %  length  of  vector  factor  for  non  random  data  vector 
n  =  0:N;  %  length  of  non-random  data  vector 

Data  Array  =  (-1) .An;  %  alternating  1,-1  square  wave  -  not  random 
Message  Length  =  length (Data  Array);  %  measures  the  length  of  the 
random  data  vector  -  not  used  with  the  non-random  data  vector 
Chipping  Matrix  =  ones  (1,64);  %  64  chips  per  bit  -  vector  of  64  ones 

%  create  the  unique  PN  code  for  transmitter  one  -  base  code  with  one 
shift 

Short  PN  =  randint ( 1 , length (Chipping  Matrix) ) *2-1 ;  %  base  short  PN 
code  with  no  shifts 

Long  PN  =  randint ( 1 , length (Chipping  Matrix) ) *2-1 ;  %  base  long  PN 
code  -  this  really  should  be  longer  than  one  bit  duration 
PN  Codel  =  circshift ( Short  PN,  [0,1])  .*Long  PN;  %  creating  the  unique 
PN  code  for  Transmitter  #ls 

%  Energy  per  bit  and  amplitude  calculations 
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Eb  =  0.5;  %  energy  per  bit  -  this  will  remain  unchanged  each  loop  - 
final  problem  should  use  0.5,  but  2  seems  to  be  working  better 
A  =  sqrt (2*Eb*Rb) ;  %  peak  amplitude 

%  this  is  the  variable  noise  loop 

SNR  =  [0,1,2,3,4,5,6,7,8,91;  %  SNR  in  DB 

for  h  =  1 : length ( SNR)  %  control  for  variable  SNR  loop 

No (h)  =  Eb/ ( 10 A ( SNR (h) / 10 ) ) ;  %  No  calculation  -  uncomment  for 
noise  loop 

%  this  is  the  transmitter/receiver  simulation  loop 
for  k  =  1 : length (Data  Array)  %  control  for  the  simulation  loop 
Chipped  Data  =  Data  Array ( k) *Chipping  Matrix;  %  repeats  the 
data  bit  Rc  times 

Mixed  Data  =  Chipped  Data.*PN  Codel;  %  mixes  the  chipped  da¬ 
ta  with  the  combined  PN  code  over  one  bit  interval 

Sampled  Data  =  upsample2 (Mixed  Data, alpha*beta) ;  %  upsamples 
the  data  to  the  sample  rate 

%  build  local  oscillator 

time  =  ( 0 :( length ( Sampled  Data) -1) ) *Ts;  %  time  vector 
Cosine  Array  =  cos (2*pi*fc*time) ;  %  cosine  vector 

%  build  the  modulated  waveform  that  will  be  transmitted 
Modulated  Signal  =  A*Sampled  Data.*Cosine  Array;  %  mixes  the 
cosine  array  with  the  sampled  data  -  result  is  the  modulated  signal 
-  prof  Kragh  changes. 

%  simulate  the  AWGN  channel  and  add  noise  to  signal 
Noise  =  sqrt (Rs*No (h) / 2 ) *randn (l,Rs) ;  %  Noise  calculation  - 
must  be  repeated  for  each  bit  in  order  to  be  independently  random 

Received_Signal  =  Modulated  Signal  +  Noise;  %  noise  is  added 
to  the  signal 

%  receive  and  process  the  signal 

Downcoverted  Signal  =  Received  Signal . *Cosine  Array;  %  mix 
with  same  cosine 

LongPN  up  =  upsample2 (Long  PN, (alpha*beta) ) ;  %  upconverts 
the  long  PN  code  to  mix  with  the  data  received  at  the  sample  rate 
Demixed  Data  =  LongPN  up . *Downcoverted  Signal;  %  mixes  the 
data  received  prior  to  LPF  input  with  the  upsampled  long  PN  code, 
b  (t) 

p  =  ones (1, Rs/Rc) *Rc/Rs;  %  averaging  factor  for  downsampling 
from  Rs  to  Rc  -  change  fc  back  to  Rc  later 

LPF_out  =  conv (Demixed  Data,p);  %  evenly  spaces  convolution 
and  averaging  -  final  LPF  output,  must  be  downsampled  still 
DFT_input  = 

LPF_out ( (0+ (alpha*beta) ) : (alpha*beta) : 64* (alpha*beta) ) ;  %  Downsam¬ 
pling  of  the  LPF  output  from  sample  rate  to  chip  rate 
G  =  fft(DFT  input);  %  DFT  of  signal  received 
C  =  fft (Short  PN) ;  %  Reference  PN  code  DFT 
CO  =  conj (C) ;  %  conjugate  of  reference  PN  code 
L  =  C0.*G;  %  mix  DFT  outputs 
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l(k, :)  =  ifft(L);  %  IDFT  of  mixer  output  -  matrix  of  data 
end  %  ends  the  transmitter/receiver  system  simulation  loop 
% - back  in  the  outer  loop 

%  make  bit  decisions  from  IDFT  output  -  idft  out  x  axis 
for  m=l : length (Data  Array) 
if  1 (m, 2 )  >=  0 

Bit  Decision (m)  =  1; 
elseif  1 (m, 2 )  <0 

Bit  Decision (m)  =  -1; 

end 

end  %  ends  the  bit  decision  loop 

%  convert  both  in  and  out  data  to  unipolar  binary 
Data  In  =  (Data  Array  +  l)/2; 

Data  Out  =  (Bit  Decision  +  l)/2; 

%  determine  error  rate 

Error  Array  =  xor (Data  Out, Data  In); 

Total  Errors  =  sum (Error  Array); 

Error  Rate (h, : )  =  Total  Errors/length (Data  In);  %  change  1  back 
to  h  after  troubleshooting 

end  %  ends  the  outer  loop  with  counter  k  -  for  variable  noise 

%  change  Error  rate  into  a  horizontal  vector  -  data  array 
Error  Rate  PN  RCVR  =  Error  Rate'; 

%export  error  rate 

save ('Error  Rate  PN  RCVR', 'Error  Rate  PN  RCVR')  %  Saves  Error  Rate 
variable  into  a  separate  file 

save ('SNR', 'SNR')  %  Saves  EbNo  array  variable  into  a  separate  file 
save('Eb', ' Eb ' )  %  saves  the  value  of  Eb  -  it  will  be  included  on  the 
plot 

%  Error  performance  will  be  plotted  in  a  separate  file. 

%  End  File 


SUPPORTING  FUNCTION  'UPSAMPLE2 '  MATLAB  CODE 

%  Name:  LCDR  Frank  Cowan 
%  Title:  upsample2  -  a  function 
%  File  Name:  upsample2.m 
%  Simulation  Project:  NPS  EE  Thesis 
%  Date  Created:  23  Apr  09 
%  MATLAB/Simulink  Version:  7.6.0,  R2008a 

%  This  function  upsamples  a  given  data  sequence  by  a  user-entered 
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factor . 

%  It  is  designed  as  a  component  of  PN  Code^RCVR  VF,  the  simulation 
of  the 

%  cyclic  PN  code  for  my  thesis. 

function  [out]  =  upsample2 (in, U) 
out  =  [ ]  ; 

for  k  =  1: length (in) 

out  =  [out  in (k) *ones (1,U) ] ; 

end 

%  End  File 
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